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Abstract. For high temperatures and densities, stellar opacities obtained from 
the Opacity Project (OP) were smaller than those obtained from the OPAL 
project. Iglesias and Rogers [Astrophys. J. 443 469 (1995)] suggested that 
the discrepancy was due to the omission by OP of important atomic inner-shell 
processes, and considered in detail results for a mixture of 6 elements: H, He, 
C, O, S and Fe. Extensive new inner-shell data have now been computed using 
the code AUTOSTRUCTURE. It is shown that the inclusion of these data in the OP 
work gives opacities for the 6-element mix which are in much closer agreement 
with those from OPAL. We also discuss a number of problems relating to the 
calculation of opacities and of equations-of-state in dense plasmas. 
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1. Introduction 

Energy is produced by nuclear reactions at the centre of a star, at a temperature 
of a few times 10 7 K, and escapes at the stellar surface. The structure of a star is 
determined by the equations for conservation of mass and of energy, an equation for 
hydrostatic equilibrium, and by the temperature gradient (see [1, 2, 3]). In regions 
for which convection does not occur, the temperature gradient is determined by the 
Rosseland-mean opacity, which is the concern of the present paper. 

In a stellar interior, an atom of a chemical clement k can exist in a number of 
ionization stages, z, and energy levels, j. Using the frequency variable u — hv / 1 (k^T) 
(where fee is the Boltzmann constant), let (Jijk(u) be the cross section for absorption§ 
or scattering of radiation by level (ijk) and let pijk be the probability of that level 
being populated. Then, the opacity cross-section for element k is 

07c (u) = SjPijkVijkju) ■ (1) 

id 

With fk the fractional abundance for element k, normalized to 

;/* = i, (2) 



E. 



For absorption processes, the correction factor for stimulated emission, [1 — exp(— u)], is included. 
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the mean cross-section per atom for a mixture of chemical elements is 

e(u) = y^JkQk{u) ■ (3) 



k 

The Ross eland- mean cross-section is ctr, where 
1 f°° 1 



F(u) du (4) 



o-r J o-{u)' 
and 

F{u) = [15/(4^ 4 )]m 4 exp(-«)/[l - cxp(-u)] 2 (5) 

(see [1, 2, 3]). Astronomers usually use opacities per unit mass. The Rosseland-mean 
opacity per unit mass is 

kr = o-r//" , (6) 

where /x is the mean atomic weight. 

The 1982 paper of Simon [4] suggested that possible errors in opacities might 
explain discrepancies between theory and observations for pulsational properties of 
stars, and provided the stimulus for two major new efforts in opacity calculations: one 
referred to as OPAL at the Lawrence Livermore National Laboratory by C. A. Iglesias 
and and F. J. Rogers; and the other referred to as The Opacity Project, OP. Early 
OPAL results [5] showed that the inclusion of very large numbers of spectrum lines 
did indeed lead to major revisions in opacities. A number of subsequent OPAL papers 
have been published, of which one of the most recent is [6] . The present paper is a 
continuation of the OP work. The calculations of atomic data were described in a 
series of papers in this journal, of which the first was [7], and OP results for opacities 
were published in [8] , to be referred to as SYMP. A collection of papers from the OP 
work, together with selected atomic data tables, have been published in book form 
[9]- 

For regions of stellar interiors with temperatures of a few times 10 5 K, which are 
of particular importance for pulsation studies, opacities from OPAL and OP are in 
close agreement and can be larger than values previously adopted, by up to factors 
of about 3. The new opacities have led to substantial improvements in the agreement 
between calculated and observed pulsation properties. However, for deeper layers, the 
results from OPAL were larger than those from OP, by amounts of 30% or so. 

Very large numbers of spectrum lines are included in the calculations, of order 
a few times 10 7 for a complete opacity run. OPAL used a parametric potential 
model [10], while OP (see [9]) used i?-matrix calculations supplemented by data from 
Kurucz [11] and from superstructure [12]. Much of the OP data includes allowance 
for configuration-interaction effects which are not included in the OPAL work. In 
their 1995 paper [13], to be referred to as IR95, Iglesias and Rogers suggested that 
the discrepancies between OPAL and OP at the higher temperatures and densities 
occurring in the deeper layers of stellar interiors could be due to the omission by 
OP of some important inner-shell transitions. We have now made new calculations, 
using the code AUTOSTRUCTURE [14], for promotions of inner-shell electrons, both via 
photoionization and photoexcitation of autoionizing states. Our results confirm the 
essential correctness of the suggestion made in IR95. The present paper describes the 
new atomic-physics calculations and presents results for the 6-element mix of IR95. 

The structure of the paper is as follows: in section 2 we discuss issues relating 
to the equation-of-state; in section 3 we detail the new inner-shell transitions that 
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we now include and the methodology that we used to describe them; in section 4 we 
briefly discuss free-free transitions: in section 5 we present and discuss our results 
for Rosseland-mean opacities; in section 6 we discuss a number of issues that arise in 
determining opacities; in section 7 we look at the opacity in the solar centre region; 
and finally in section 8, we give a brief summary. 

2. The equation-of-state (EOS) 

The fraction of element k in ionization stage i is 

4>%k = ^Pijk ■ (7) 

3 

We may put 

Pijk = 4>ik x gtjkWijk exp[-E ijk / (k B T)]/U ik (8) 

where, for level {ijk), gijk is statistical weight, Eijk the total energy and W^k an 
occupation probability. The internal partition function is 

Uik = ^gijkWi-jk exp[-E ijk /(k B T)] . (9) 

3 

If all Wijk were set to unity, the summation for Uik would be divergent. Due to 
interactions with particles and fields of the surrounding plasma, states of sufficiently 
high energy have only small probabilities of being occupied and, hence, small values 
of Wij k . 

2.1. Occupation probabilities 

The methods used in the OP work for the calculation of the Wijk are described by 
Hummer and Mihalas in [15]. Let F be the ion micro-field and P(F) the micro-field 
distribution (MFD): P(F)dF is the probability of F being in the range dF, with 
/ °° P(F) dF = 1. Hummer and Mihalas define critical fields, Fijk, such that species 
(ijk) can only exist in an environment with F < F^k, giving 

W ijk = ( Fl3k P{F)dF. (10) 
Jo 

The use of equation (10) gives W ~ 1 for low densities but W to be small for high states 
and high densities. In [15] Hummer and Mihalas used the Holtsmark MFD, which does 
not make full allowance for correlations between particles nor for the charge on the ion 
(ijk). They also introduced a simple analytical approximation (equation (4.70) of [15]) 
which gave results in close agreement with those from (10) at low densities but a more 
rapid exponential decrease at high densities. A similar exponential form is obtained 
using a 'hard-sphere' approximation (see section III (a) of [15]). The analytical form, 
(4.70) of [15], was used in SYMP. 

Iglesias and Rogers, in IR95, noted that the Holtsmark distribution would not 
be a good approximation at the higher densities of stellar interiors and Nayfonov et 
al. [16] subsequently obtained improved expressions for the Wijk using the micro-field 
distributions of Hooper [17], which are also used in the OP calculations of line-profiles 
for hydrogenic ions [18]. Following Nayfonov et al, we refer to the equation of state 
using their MFD as the Q-EOS. Their Wijk are fairly close to those obtained using 
the APEX distribution of IR95 and, with one modification to be described below, will 
be used in the present work. 
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2.2. The ionization equilibrium 



For simplicity of presentation, we here omit specification of the clement index k and 
neglect all refinements such as allowance for electron degeneracy. The conclusions 
reached remain valid when all necessary refinements are included. We take the 
ionization index i to be equal to the number of bound electrons: i = for the bare 
nucleus and i = Z for a neutral atom with nuclear charge Z . 

The ratio of ionization fractions in successive stages is (see [7] or [15]) 



i x U c 



3/2 



4>i-\ u, 

where N c is the electron density, 

' m c k B T \ 
2irh 2 J 

and m e is the mass of an electron. 

Since Uq = 1, the fraction in stage i relative to the fraction of bare nuclei is 



(11) 



(12) 



= U t x 



'Ne 



(13) 



2.3. Pressure ionization 

Now consider fixed T and increasing N c . Initially, <pi / (f>o will increase with increasing 
N c due to the factor [N c /U c } 1 in (13). That effect is pressure recombination. 
Eventually, for large densities, the Wij in Ui become small and, if they become 
sufficiently small, pressure recombination can be followed by pressure ionization. 

For W calculated using (10), in the limit of N e large one obtains W oc N~ 2 using 

3/2 

the Holtsmark MFD and W oc N e ' with the Q form. In either case it is seen that 
pressure ionization does not occur for i > 2. That result would appear to be quite 
unphysical, since it would imply that all atoms are in states which have very small 
occupation probabilities. Pressure ionization does occur, for all i, if one uses a 'hard- 
sphere' model for W or, as in SYMP, the approximation of equation (4.70) of [15], 
giving exponential decreases of W with increasing N e . 

We conclude that equation (10) can be expected to give reasonable values for W 
when W is not very small, but that it does not give a sufficiently rapid decrease of W 
in the limit of high densities. We adopt the expedient of introducing a critical value 
W c of W; use the value of W from equation (10) if it is greater than W c ; and take 
W = if equation (10) gives W < W c . For quite a wide range of values of W c , results 
for opacities are found to be insensitive to the value of W c adopted. The final results 
reported in the present paper are obtained with W c — 10 -3 . 

The need for introducing a cut-off in W is illustrated in Figure 1 for carbon 
at rather high temperature (log(T) = 7.5) and very high densities (log(7V ) = 27.0, 
27.5 and 28.0). Figure 1 (a) shows ground-state occupation probabilities W against 
ionization stage i. Neutral carbon, i = 6, is seen to have very small values of W which 
decrease with increasing N e . Figure 1 (b) shows ionization fractions calculated without 
a cut-off in W. There is seen to be an abrupt change from the case of log(iV c ) = 27.0 
with carbon nearly fully-ionized (ionization fraction 0.98 for i = 0) to the case of 
\og(N c ) — 28.0 with carbon nearly fully-neutral (ionization fraction 0.95 for i = 6). 
The result for log(iV ) = 28 is clearly nonsensical, with nearly all of the carbon in a 
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state with occupation probability of W ~ 1CT 12 ! With a cut-off of W c = 1CT 3 we 
obtain a result which is much more plausible: for log(iV c ) = 27.0, ionization fractions 
very similar to those for W c = for i < 3 and equal to zero for i > 3; and for 
log(iV e ) = 27.5 and 28.0 the carbon is fully ionized. 

2.4. The W n of OPAL 

Bethc and Salpeter [19] give expressions for (r 3 ; ) for hydrogenic ions of charge Z. We 
put 

<^> = (4)x>+i)<^> ( 14 ) 
vn * 1=0 

to obtain 

(r 3 ) = [n 2 /(8Z 3 )] [21n 4 + 35n 2 + 4] . (15) 
The mean volume for state n can be defined as 

K-(47r/3)(r 3 ). (16) 

Let N e be the electron density, N a be the atom density and p be the mass density. 
From Table 1 and Figure 1 of IR95, N e = 5.0 x 10 21 cur 3 and, for log(p) = -2 and 
the 6-element mix (see section 5.2), N a = 4.7 x 10 21 cm~ 3 . The total particle density 
is then N = N c + N a = 9.7 x 10 21 cm~ 3 . The average volume occupied by a particle 
(electron or atom) is 1/N. 

Table 1 gives, for hydrogenic carbon and the case of Table 1 of IR95 (log(T) = 6, 
log(p) = —2): values of NV n ; VF„(OPAL) from IR95; and W„(Q) from formulae given 
in [16]. The values of W n (Q) are fairly close to the values of W n (APEX), as given in 
IR95. 

In both OP and OPAL, optical properties (oscillator strengths and photoion- 
ization cross sections) are calculated for unperturbed atomic states, and W n is the 
probability of state n being unperturbed. With NV n > 1, one would expect to find at 
least one other plasma particle within the volume V n , giving a state which is markedly 
perturbed and which one would therefore expect to have a small value of W n . By that 
criterion, some values of W„(OPAL) seem to be surprisingly large, particularly those 
for n = 4, 5 and 6. Values of W n (Q) appear to be more reasonable. 



3. Inner-shell atomic physics 

The original OP work utilized the i?-matrix method which uses wavefunction 
expansions of the type 

* = ^^V„0„, (17) 

n 

where the ip n are functions for atomic 'target' states, the 8 n are functions for an added 
electron, and A is an anti-symmetrization operator. In that method, photoionization 
and autoionization are treated as a single quantum-mechanical process. Figure 2 gives, 
as an example, the cross section for transitions from the ls2s 1 S state of Fe 24+ . Below 
the threshold for ejection of the Is electron, the cross section shows autoionization 
features due to processes 

ls2s + hv ->• 2snp ->■ Is + e" . (18) 
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However, use of the i?-matrix method would not be practicable for the computation 
of the atomic data required for the present work, for two reasons: (a) in many 
cases, the number of channels n required in equation (17) would be prohibitive (our 
experience with the RmaX work on inner-shell X-ray processes [20] shows that a few 
Li- and Be- like ions are the most that could be treated with a reasonable timescale); 
(b) it would be difficult to allow for pressure-broadening of the autoionization 
features. We therefore use a perturbative approach, as implemented in the program 

AUTOSTRUCTURE [14, 24]. 

We note that Figure 2 shows, qualitatively, the effects of interference 
between autoionization features and the background continuum. However, detailed 
quantitative studies [21] show that, on averaging over resonance profiles, this 
interference is a very small effect and it can safely be neglected for our purposes. 
This is the independent processes approximation. In the time-reversed case, this 
corresponds to treating dielectronic and radiative recombination separately. The 
second approximation required by our perturbative approach is an isolated resonance 
treatment of the autoionizing features. The effect of interacting resonances has been 
investigated [21] for the reverse process of dielectronic recombination, and it also can 
safely be neglected for our purposes. 



3.1. Photoexcitation 

The downward probability rate for a radiative transition from an upper state u to a 
lower state I is given by 

where g u is the statistical weight for the upper level, c is the speed of light, u) = 2kv 
where v is the photon frequency, and S\ u is the bound-bound line-strength, as 
defined in [22]. Eissner et al [23] give expressions for Si u in multi-configuration 
LS- and intermediate-coupling (their equations (115) and (117), respectively): their 
code, superstructure, gives numerical values for Si u in atomic units. The code 
AUTOSTRUCTURE, which incorporates superstructure, gives the following data: 
(a) the line-centre frequency; (b) the rates {g u / 9l)-^ T u ^i', ( c ) A T U = J2i ^-u^n the 
total probability rate for radiative decay of the upper level u; (d) A*^ m and = 
J2m ^u^rm the partial and total autoionization decay probability rates, respectively. 
Multi-configuration LS- and intermediate-coupling expressions for A a a _ ¥m are given in 
[14] (equations (2.2) and (2.4), respectively). The rates are all in s _1 . 
The oscillator-strength for the I — > u transition is 

= m A '-' m 

and the absorption (photoexcitation) cross section for an autoionizing feature is 

oH u = — (21) 



m e c 

where <p(w) is the line-profile factor. On neglecting thermal Doppler broadening, the 
profile is 

<p(u>) - (T/2tt)/[(w - ^o) 2 + (r/2) 2 ] , (22) 
where loq is the line-centre angular frequency and 

r=(Al + A* u )+AZ + Tl, (23) 
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where expressions for the pressure-broadening contribution, rj^, are given in [35] (a 
different form is used for hydrogenic lines, see [18]). The profile (22) is then convolved 
with that for Doppler broadening to give a Voigt profile. 

For more detailed (collisional-radiative) modelling purposes we require to follow 
the break-up of an autoionizing state. For example, the photoexcitation cross section 
in equation (21) is multiplied by branching ratios for radiative decay, A T u ^ l /(A I u +A^ l ), 
and by Auger yields for autoionizing decay, A^ m / (A T U + A*). 

3.2. Photoionization 

The direct photoionization cross-section from an initial state I of an (N + l)-electron 
atom to a final state u of an TV-electron ion plus ejected electron is given by [7] 

<^ = 7^£W (24) 
gi 6c 

where gi is the statistical weight of the initial level and S u i is the bound-free line- 
strength with the final-state continuum wavefunction normalized per unit energy. The 
code AUTOSTRUCTURE uses a different continuum normalization and gives numerical 
values for the bound-free line strength which are equal to it/ 2 times S u i in atomic units. 
The same expressions for S u i that were used for bound-bound transitions are used 
for bound-free transitions, except that the final 'active' bound-state wavefunction is 
replaced by a continuum distorted-wave function which does not contain any resonance 
structure. This use of distorted- waves is a good approximation for atoms that are a few 
times ionized and is the final approximation employed by our perturbative approach. 

3.3. AUTOSTRUCTURE, some details 

3.3.1. Angular algebra AUTOSTRUCTURE [14, 24] incorporates SUPERSTRUCTURE 
[23] and the angular algebra required to calculate the preceding atomic data is no 
more than that which is generated by superstructure to determine energy levels 
and radiative rates. (The angular algebra required for autoionization rates is that 
which follows from the H operator which determines the structure.) However, for the 
complex inner-shell processes considered here, which can give rise to configurations 
with thousands of terms, we found it necessary to re-write the angular algebra 
code. Specifically, because of historical memory limitations, the algebraic Slater- 
state interactions were determined between each LS term or J level [23]. However, 
the Slater-state interaction depends only trivially on the configuration and can 
be generated much more efficiently between symmetry groups. With complex 
configurations there is a high degree of algebraic term and level degeneracy, both 
within and between configurations. Re-coupling by LS or LSJ symmetry groups 
reduces the overall time spent on the largest scale jobs by a factor of 30-40 (LS) or 
100-200 (LSJ). 

3.3.2. The Hamiltonian matrix superstructure determines multi-configuration 
eigenenergies and eigenvectors by diagonalization of the Hamiltonian matrix for each 
SLw or Jir symmetry group. AUTOSTRUCTURE further partitions the problem by 
(N + l)-electron bound and iV-electron (plus continuum) configurations — the N- 
and (N + l)-clectron Hamiltonians arc diagonalizcd separately. The bound-free 
Hamiltonian interaction between the two is treated as a perturbation — this leads 
simply to the autoionization rate [14] . 
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The mass-velocity and Darwin operators of the Breit-Pauli Hamiltonian [23] are 
also added to the usual non-relativistic Hamiltonian for the determination of our LS- 
coupling atomic structure. This results-in transition energies between terms which are 
in good agreement with those obtained from using intermediate coupling, on averaging- 
over fine-structure. 

3.4- Application to the 6-element mix 

The JT-shell processes required are of the form: 

ls q 2l p n'l' + hv^ ls q - 1 2Pn'l' + e" (25) 
for photoionization and 

ls q 2l p n'l' -\- hv — ► ls q - 1 2l p n'l'n"l" -» la q - 1 2l p n"T' + e" (26) 

for photoexcitation-autoionization, where 2l p stands for 2s s 2p* with p = s + t. 
Calculations are made for q = 1, 2 for p = and q = 2 for p > 0. The values 
of p depend on the element, for example, up to p = 7 for iron but only p = 1 for 
carbon. This depends on the importance of its contribution to the opacity. We use 
n',n",n"' — 2 to 6, for all allowed I', I", I'". The contributions from higher-n in (26) 
are obtained by matching onto the results of (25). 
The L-shell processes required are 

2l q 3l' p n"l" + hu^ 2l q - 1 3l' p n"l" + e" (27) 

and 

2l q 3l' p n"l" + hu^ 2l q -Hl' p n"l"n"'l"' -» 2l q - 1 3l' p n iv l iv + e~ , (28) 

where 3l' p stands for for 3s s 3p*3d" withp = s + t + u. Calculations are made for q = 1 
to 8 for p = and q = 8 for p > 0; e.g., up to p = 2 for iron. 
The M-shell processes required are simpler: 

3l q n'l' + hv -> 3Z g -Vr + e" (29) 

and 

3l q n'l' + hv^ 3l q - 1 n'l'n"l" -> ?A q - 1 n"'l l " + e" , (30) 

and are included only for iron, with q = 1 to 3. 

To attain as much accuracy as possible for consistency with the existing i?-matrix 
data, our perturbative calculations retain configuration interaction within the ./V- 
clcctron (core) complex and the (N + l)-electron complex (for n' — n"): for example, 
we retain interactions between states such as 2s 2 1 S and 2p 2 X S. For most cases, 
calculations are made both in LS'-coupling and in intermediate coupling (including 
Breit-Pauli terms). Even when not significant for inner-shell contributions to opacity, 
this resolution is required for modelling non-LTE photoionized plasmas. 

In table 2, we summarize the inner-shell calculations that we have carried-out. 
We have examined the convergence with respect to the inclusion of inner-shell data 
and estimate that the inclusion of further data would change the Rosseland means by 
less than 1%. 

The new data resulted-in an additional 2 575 458 level-resolved photoexcitation 
lines and 187351 total photoionization cross sections, i.e. summed-over all final states, 
contributing to the OP opacities. Much more data are archived, viz. final-state 
resolved photoionization cross sections and data for transitions and lines than were 
already included in the original OP work. 
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3.5. Database issues 

Given that we have already generated, and will be generating, large amounts of 
atomic data, some thought has been given to how it might be effectively archived 
for applications other than the one at hand. In particular, the Atomic Data and 
Analysis Structure (ADAS) [25] has long had the capability of handling radiation 
fields [26] but has, until recently, concentrated on utilizing escape factors [27]. The 
collisional-radiativc modelling of finite-density non-LTE photoionized plasmas with 
ADAS requires that we archive final-state resolved photoionization data, summed- 
over final channel angular momenta. We have written codes (adaspe and ADASPi) 
to process the large amount of energy levels, radiative rates, autoionization rates and 
photoionization cross sections produced by AUTOSTRUCTURE and have defined such a 
suitable archive data structure of final-state resolved photoexcitation-autoionization 
and direct photoionization data, specifically, adf38 and adf39 ADAS data formats 
[25] . This partial data is then further reduced to total photoionization data for use by 
opacity calculations. For completeness, outer-shell photoionization data is archived in 
the adf39 files as well as inner-shell, but it is not used in the work reported-on here. 

4. Free free transitions 

Contributions from free-free transitions are, in most cases, calculated in a hydrogenic 
approximation using the code of [28] which employs fits to accurate results from [29] . 
Allowance for Debye screening can be of importance at high density and is included if 
it gives increases in Rosseland-means by more than 0.1%. The screening contribution 
is calculated using the Born-Elwert theory ([19, 30, 31]). It is checked that, without 
screening, the Born-Elwert approximate gives agreement with the accurate results 
from [28] and [29] to better than 1% for all cases considered in the present work. For 
log(i?) < — 1 and log(T) < 7.7 it is found that inclusion of screening for free-free 
transitions never changes the Rosseland-mean by more than 2 or 3 %. 

The correction for Debye screening was not included in the work of SYMP but is 
included in the present work. 

5. Results for Rosseland-mean opacities 

5.1. Use of the variable log(i?) 

It is convenient to use the variable 

R = p/Tl , (31) 

where p is the mass density in g cm -3 and T 6 = 10~ 6 x T, with T in K. For a solar 
mix, comparisons of log(«R,) from OP and OPAL for log(i£) = — 1 to — 6 were given in 
Figure 15 of SYMP, which showed the OPAL opacities to be larger than those from 
OP for larger values of log(ii) with log(T) > 5.5. 

5.2. The 6-element mix 

Iglesias and Rogers in IR95 considered the importance of inner-shell transitions for a 
mixture of 6 elements (H, He, C, O, S and Fc) with abundances chosen to be such as 
to give opacities similar to those for the complete solar mix. The adopted number- 
fractions, ft, are given in Table 3. 
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Paper IR95 gave results for one temperature-density point, log(T) = 6 and 
log(i?) = —2 giving log(p) = —2. OPAL calculations were made both with and 
without the inclusion of inner-shell processes. The inclusion of those processes was 
found to give an increase in kr by 30%. 

5.3. OP results both with and without inner-shell transitions 

Figure 3 gives OP values of log(fCR) for the 6-element mix, both with the inclusion of 
the inner-shell data discussed in Section 3 and without. The results without inner- 
shell data are essentially the same as those given in SYMP, but include some fairly 
minor improvements and use of the Q-EOS. Figure 4 shows <51og(«R), the increase 
in log(KR) due to inclusion of the inner-shell data. For log(T) = 6, log(i?) = —2 wc 
obtain an increase in kr by 31%, in close agreement with the increase obtained in 
IR95. 

Changes in kr due to changes in the EOS are much smaller than those due to 
the inclusion of inner-shell data: thus for log(T) = 6, log(i?) = —2, use of the Q-EOS 
in place of the EOS used in SYMP reduces kr by only 1.6% (in both cases, without 
inner-shell data). 

54. Results from OP and OPAL 

Figure 5 shows values of log(KR) for the 6-element mix from OPAL (data obtained 
from the OPAL website [32]) and from the present OP work with the inclusion of 
inner-shell data. It is seen that, for the larger values of log(i?) where the inner-shell 
data is important, OP is now in close agreement with OPAL. 

There are some remaining differences between OP and OPAL at smaller values of 
log(T) and log(i?) where inner-shell data are not important. A feature at log(T) ~ 5.2, 
often referred to as the 'Z-bump', was first identified by Iglesias, Rogers and Wilson 
[33] as due to inclusion of very large numbers of M-shell transitions for various 
ionization stages of iron. For log(-R) < —3 there are seen to be some differences 
between OPAL and OP in the vicinity of that feature. Figure 6 shows, on a much 
more expanded scale, percentage difference for kr, (OP— OPAL), for log(i?) = —3 and 
—6. For log(-R) = —3 the Z-bump feature from OP is seen to be broader than that 
from OPAL, and to have a peak value lower by about 10%. For log(i?) = —6 the OP 
peak value is lower by 20%. We plan to consider these differences further in a later 
paper. 

6. Discussion 

IR95 was mainly concerned with illustrative results for the case of the 6-element mix 
and log(T) = 6, log(p) = —2. Plots of monochromatic opacities were given for C, S, 
Fe and for the mixture. Similar plots were obtained in the course of the present work. 

6.1. Density dependence 

It is seen from Figure 4 that the inner-shell contributions are most important at 
high densities. That can be understood by considering the case of high temperatures 
(similar considerations apply for other regions). At high temperatures, say log(T) > 
7.0, the most important transitions are for iron K-sheW. At the very lowest densities 
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considered, the iron is almost fully ionized and there are no if-shell contributions. 
Then, as the density increases, one has if-shell transitions of the type 

ls + hv->np, (32) 

and 

Is 2 + hv — > ls«;p , (33) 

where k = n for transitions to a bound state and n = k for photoionization. Such 
transitions are included in the original OP calculations. With a further increase in 
density, one has states of the type ls 2 C (where 'C stands for states of outer electrons) 
and if-shell transitions of the type 

ls 2 C + ^^ lsCVtp. (34) 

Such transitions were not included in the original OP work but arc included in the 
present work (see section 3). 

6.2. Line profiles 

Pressure broadening of all spectrum lines should be included, for both transitions to 
true bound states and for those to autoionizing states. 

The importance of pressure-broadening for the iron if -shell is illustrated in Figure 
7, which gives plots of a(u) for iron for log(T) = 7 and log(iV c ) = 22.5 and 25: for 
the 6-element mix, those values of N c correspond to \og(R) = —4.21 and —1.71. It is 
seen that, for the low-density case, there are very many resolved spectrum lines but 
that at higher densities pressure broadening leads to the lines being almost completely 
blended. We use the pressure-broadening theory from [35] and [18] and line-blending 
theory from [36]. 

We note that, beyond the K-edge, the cross sections are nearly the same for both 
densities, since both depend mainly on if-shell ionization. At the lower energies, say 
u < 7, the background opacity is mainly due to free-free transitions, giving cross- 
sections per atom proportional to the electron density. 

6.3. Autoionization 

The AUTOSTRUCTURE calculations give line widths due to radiation damping and to 
autoionization. We find, in practice, that inclusion of those contributions to the line 
profiles is never very important since widths due to pressure broadening are generally 
much larger: omission of autoionization widths would never give errors in Rosseland- 
means larger than one on two percent. 

6.4- Fine- structure 

In both the OPAL and OP work, it was found that inclusion of fine-structure for 
outer-shell transitions could be of importance for the calculation of Rosscland-mcans 
(see [34] and SYMP). Most of our AUTOSTRUCTURE calculations were made both in 
LS'-coupling (no fine-structure) and in intermediate coupling (with fine-structure). 
Test runs showed that inclusion of fine-structure for the inner-shell transitions never 
increased Rosseland-means by as much as 1%, which was, again, a consequence of the 
importance of pressure-broadening. All final results as reported in the present paper 
were made with omission of inner-shell fine-structure. 
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7. The solar centre region 

Convection occurs throughout much of the solar interior and in those regions a precise 
knowledge of the Rosseland-mean opacity is not of great importance. However, 
convection does not occur in the deepest layers of the solar interior. Knowledge of the 
Rosseland-mean for the centre region is of importance for the construction of solar 
models, which can be tested against data from helioseismology. 

At the solar centre, models give log(T) = 7.196, log(p) — 2.179 and log(i?) = 
— 1.409 (see [37, 38]). Opacities calculated using a number of different codes [39] show 
quite a large scatter. The most accurate values currently available are undoubtedly 
those from OPAL [40] . For the 6-clcment mix, opacities for the centre region from OP 
are a little larger than those from OPAL: by 2.3% at log(T) = 7.2 and log(i?) = —1.5. 

The present work permits some further discussion of the solar-centre problem. 
In the OP work we use a mesh of values of log(T) and log(iV ). We select a mesh- 
point log(T) = 7.2, log(N c ) = 26 giving, for the 6-element mix, log( i o) = 2.294 and 
log(i?) = —1.306. Figure 8 shows log(a) for the mixture at that point. There are 
three main contributions to the centre opacity. 

7.1. Free-free 

The free-free contribution is mainly due to electron collisions with H + and Hc +2 . The 
free-free cross section behaves like u~ 2 and the process therefore dominates at low 
frequencies. We include Debye screening (see Section 4) which modifies the solar- 
centre opacities by about 1 or 2 per cent. OPAL include some further refinements (see 
[13]) but they are not likely to be of much importance. 

7.2. Electron scattering 

At lower densities, the electron-scattering cross section is equal to the Thomson cross 
section, and is independent of frequency. For higher densities, the cross section is 
modified by plasma collective effects. Both OP and OPAL use the theory of Boercker 
[41], and also allow for relativistic corrections. 

7.3. Atomic transitions 

In Figure 8 there is a just-discernible feature at u ~ 2 due to sulphur if-shell 
transitions, and a much more conspicuous one at u ~ 5 due to iron if-shell transitions. 
For hydrogenic iron, the Ly Q line is at u — 5.10 and the Lyman continuum starts at 
u = 6.80. Table 3 gives ionization fractions and ground-state occupation probabilities 
for iron: all stages up to i = 9 are seen to contribute to the ^-feature. It is seen from 
Figure 7 that the high K-sheW lines are completely blended. It may be noted that the 
mean opacity will not be sensitive to the exact distribution amongst ionization stages 
(the 4>i of Table 4) since the cross sections for promotion of Is 2 electrons will be much 
the same for the different stages. 

7.4- Use of a different frequency variable 
Equation (4) may be replaced by 




(35) 
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where 

v(u) = [ F(u)du (36) 

Ju=0 

and j) max = v(u — > oo) — numerical integrations give w max = 1.0553. Figure 9 shows 
l/cr plotted against v for the case of Figure 7. The advantage of using Figure 9 is that 
it shows the sensitivity of 1/<tr, to the various features in o~(u). There are seen to be 
rather small contributions from regions of small u where a(u) is large but much more 
important contributions from large u where a(u) is small. 

8. Summary 

Rosseland-mean opacities from the Opacity Project, OP, were originally found to 
be smaller than those from the OPAL project at high temperatures and high densities. 
Iglesias and Rogers, in IR95 [13], discussed the case of log(T) = 6, log(p) = —2 where 
kr(OPAL) was larger than kr(OP) by about 30%. They made two criticisms of the 
OP work: (a) OP calculated occupation probabilities W using a Holtsmark MFD 
where, at high densities, it is a poor approximation; (b) OP omitted some important 
inner-shell atomic data. 

8.1. Occupation probabilities 

In the present work we calculated values of W using expressions from Nayfonov et al 
[16], who used an MFD theory that is valid to high densities. It is shown that the 
expressions for W given in [15] and [16] are not valid in the limit of very high densities: 
a simple expedient removes that difficulty. It is found that final results for Rosseland- 
mean opacities are not very sensitive to the adopted occupation probabilities. 

It is noted that the values of W obtained from the OPAL work seem to be rather 
surprisingly large for more highly-excited states. 

8.2. Inner-shell data 

New opacity calculations have been made for the 6-element mix introduced in IR95. 
Extensive inner-shell atomic data were computed using the code AUTOSTRUCTURE for 
the 6-element mix introduced in IR95. It is shown that inclusion of those data removes 
all major diffcrcncics between the OP and OPAL work. 

8.3. Future work 

In the work described in SYMP, opacities were calculated for 17 cosmically- abundant 
chemical elements. Work is now in progress to obtain inner-shell data for those 
elements not included in the present work. 
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Table 1. Occupation probabilities (W n ) and 

average volumes (V n ) for C 5 + at 

T = 10 6 K and N = 9.7 X 10 21 cm- 3 . 



n 


NV„ 




VMOPAL) 


W„(Q) 


1 


2.08(- 


-4) a 


1.000 


1.000 


2 


6.67(- 


-3) 


0.996 


0.997 


3 


6.31(- 


-2) 


0.995 


0.967 


4 


0.330 




0.995 


0.705 


5 


1.216 




0.914 


0.154 


6 


3.562 




0.527 


1.58(-2) 


7 


8.875 




0.162 


1.87(-3) 


8 


1.96(4 


-1) 


2.37(-2) 


2.82(-4) 


9 


3.96^ 


-1) 


2.23(-3) 


5.22(-5) 



a 2.08(-4) = 2.08 x 10" 4 . 



Table 2. Inner-shell transitions considered: 

Z is the nuclear charge of the clement and N the 

number of electrons left on the final ion. 



z 


N 


Shells 


Z 


N 


Shells 




2 


1 


K 










6 


1 


K 


26 


1 


K 






2 


K 




2 


K 






3 


K L 




3 


K L 












4 


K L 




8 


1 


K 




5 


K L 






2 


K 




6 


K L 






3 


K L 




7 


K L 












8 


K L 




16 


1 


K 




9 


K L 






2 


K 




10 


L 






3 


K L 




11 


L 


M 




4 


K L 




12 


L 


M 




5 


L 




13 




M 




6 


L 












7 


L 












8 


L 












9 


L 












10 


L 
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Table 3. Number fractions 
for the 6-element (k) mix of 
Iglcsias and Rogers [13]. 



k 


fk 




H 


9.071(- 


-1) 


He 


9.137(- 


-2) 


C 


4.859(- 


-4) 


O 


9.503(- 


-4) 


s 


9.526(- 


-5) 


Fe 


3.632(- 


-5) 



Table 4. Ionization fractions (<f>i) and 
ground-state occupation probabilities (Wo) 
for iron for log(T) = 7.2 and log(7V c ) = 26. 



i 


4>i 


Wo 





0.000 


1.000 


1 


0.001 


1.000 


2 


0.039 


1.000 


3 


0.155 


0.959 


4 


0.289 


0.954 


5 


0.282 


0.941 


6 


0.162 


0.930 


7 


0.058 


0.918 


8 


0.012 


0.894 


9 


0.002 


0.872 


10 


0.000 


0.847 
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Figure 1. Carbon at log(T) = 7.5 and log(JV e ) = 27.0, 27.5 and 28.0. (a) Ground- 
state occupation probabilities, W, against ionization stage i. (b) Ionization 
fractions, (pi, calculated without a cut-off, W c , in W. 




Figure 2. The atomic cross-section, Cjjfc, for photoionization from Fe 24+ ls2s 1 S, 
in the vicinity of the if-edge. Cross section in atomic units. 



Inner-shell contributions to opacity 



19 





Figure 4. 51og(Kpt), the change in Rosseland-mean opacity which results from 
the inclusion of inner-shell data. 
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LOG(T) 



Figure 5. Rosseland-mean opacities for the 6-element mix, with the inclusion 
of inner-shell data: full lines, OP, present work; dashed lines, OPAL, from [32]. 
Opacities in cgs units. 
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Figure 7. The opacity cross-section, a, for iron in the vicinity of the if-edge, at 
log(T) = 7: cross section in atomic units. Pressure broadening included, 
(a) log(Ar o ) = 22.5, and (b) log(AT e ) = 25. 
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Figure 9. 1/cr against v, denned by equation (35), for the case of Figure 



